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Abstract: The phase diagram of quark gluon plasma (QGP) formed at a very 
early stage just after the heavy ion collision is obtained by using a holographic dual 
model for the heavy ion collision. In this dual model colliding ions are described by 
the charged shock gravitational waves. Points on the phase diagram correspond to 
the QGP or hadronic matter with given temperatures and chemical potentials. The 
phase of QGP in dual terms is related to the case when the collision of shock waves 
leads to formation of trapped surface. Hadronic matter and other confined states 
correspond to the absence of trapped surface after collision. 

A multiplicity of the ion collision process is estimated in the dual language as an 
area of the trapped surface. We show that a non-zero chemical potential reduces the 
multiplicity. To plot the phase diagram we use two different dual models of colliding 
ions, the pointlike and the wall shock waves, and find qualitative agreement of the 
results. 
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1. Introduction 



For the last decade, since the publication of fascinating papers |], @, [|, it was realized 
that supersymmetric and non-supersymmetric theories in the strong coupling limit in 
principle could be pretty close in their properties Q. The AdS/CFT correspondence, 
which appeared as a formal duality between the M = 4 super Yang- Mills theory and a 
quantum gravity in AdS background, has become powerful tool for studying various 
properties of real physical systems in the strong coupling ||. 

Important branch of these investigations is the analysis of the Quark Gluon 
Plasma (QGP) from the point of view of ArfS'-holography, see for example, review 
|J. These applications of the AdS/CFT correspondence to strongly coupled QGP 
have been mostly related to equilibrium properties of the plasma, or to its' kinet- 
ics/hydrodynamics near the equilibrium. 

A particular application of AdS/CFT to the strongly coupled QGP, is the anal- 
ysis of thermalization of matter and early entropy production instantly after the 
collision of relativistic heavy ions. RHIC experiments have shown that a QGP forms 
at a very early stage just after the heavy ion collision, i.e. a rapid thermalization 
occurs, and QGP produced in RHIC is believed to be strongly coupled as evidenced 
by its rapid equilibration. Strong collective flows well reproduced by hydrodynam- 
ics, and strong jet quenching f?|, §, ||. This obviously requires a calculation of the 
strongly coupled field theory in non-equilibrium process. 

Not long ago Gubser, Yarom and Pufu [[Uj have proposed the gravitational shock 
wave in AdSs as a possible holographic dual for the heavy ion and have related the 
area of the trapped surface formed in a collision of such waves to the entropy of 
matter formed after collision of heavy ions. Early papers where has been mentioned 
an analogy between colliding heavy ions and colliding gravitational shock waves in 
anti-de Sitter space include [Q-[]l6). This AdS-holographic model has been also used 
to find the stress-energy tensor of the QGP formed by ion collision. In accordance 
with AdS / CFT dictonary this stress-energy tensor is dual to the metric of spacetime 



after collision of shock waves [16 



The main result of [11, [7], confirmed by numerical calculations performed in 



1~9| , is that in the limit of a very large collision energy E the multiplicity (the 



entropy S) grows as 

S > CE 2/3 , (1.1) 

C is a numerical factor (see Sect. 2. 1.4 ). 

Alvarez-Gaume, Gomez, Sabeo Vera, Tavanfar, and Vazquez-Mozoand |20| have 
considered central collision of shock waves sourced by a nontrivial matter distribution 
in the transverse space and they have found critical phenomenon occurring as the 
shock wave reaches some diluteness limit. This criticality may be related to criticality 
found in ]n|. The numerical results of |]1J show the existence of a simple scaling 
relation between the critical impact parameter and the energy of colliding waves. 
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The size of colliding nuclei is introduced via the distance of those objects from 
the boundary along the holographic coordinate z. 

The model of infinite homogenous wall has been proposed and analyzed by 



Shuryak and Lin [18|. The advantage of this model is the essential simplicity of 
calculations. However, the legitimacy of these calculations requires some additional 
examinations (see our discussion in Sect. |2.2| ). 

In heavy ion collisions not only the energy per nucleus is important variable. One 
can try to associate different nuclei with different kinds of shock waves. There are 
several proposal in literature on this subject. For example, in [[H] the holographic 
model with cutting off the UV part of the bulk geometry, has been proposed. For- 
mation of trapped surfaces (TS) in head-on collisions of charged shock waves in the 



(A)dS background has been considered in E2j and it has been shown that the forma- 



tion of trapped surfaces on the past light cone is only possible when charge is below 
certain critical value - situation similar to the collision of two ultrarelativistic charges 
in Minkowski space-time |]23|| . This critical value depends on the energy of colliding 
particles and the value of a cosmo logical constant. The formation of trapped surfaces 
in head-on collisions of shock waves in gravitational theories with more complicated 
bulk dynamics, in particular with the Einstein-dilaton dynamics, pretended to de- 



scribe holographic physics that is closer to QCD than the pure AdS theory [p4| , ^1 
has been considered recently by Kiritsis and Taliotis pSJ 1 and they have found that 
the multiplicity grows as 

S>E°- 2 \ (1.2) 

that is rather close to the experimental data. 

In this paper we propose to incorporate the study of collisions of charged shock 
gravitational waves [22[] into the description of colliding nuclei with non-zero baryon 



chemical potential. In the holographic context, the chemical potential of strongly 
coupled QGP on the gravity side is related to temporal component A t of the U(l) 



gauge field |27|]-[[34j]. The asymptotic value of this gauge field component in the bulk 



is interpreted as the chemical potential in the gauge theory 

H = A t \ boundary • 

(1.3) 

We use the same identification ( |1.3|) for colliding ions. It would be interesting to 
perform calculations for the off-center collision of charged gravitational waves or 
generally smeared charged shock waves. Postponing this problem for further inves- 
tigations, here we consider the head-on collision of charged point shock waves and 
charged wall shock waves. This will give us the holographic picture for QPG phase 
diagram in the first moment after collisions of heavy ions. This phase diagrams, 
chemical potential (charge) fi versus temperature (energy) T, are displayed in Fig. 5 
and Fig. 11. The colored lines separate the TS-phase from the phase free of TS. Let 



1 Collision of dilatonic shock waves in the fiat background has been considered in 
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us note that the obtained diagrams differ from the phase diagram for equilibrium 
QGP (see Fig.|in Sect. |T|). 

The paper is organized as follows. In Sect. 2 we present our set up of the problem. 
In Sect. 2. 1.1 we describe the role of black holes in AdS/CFT description of strongly 
coupled QGP. In Sect. 2.1.2. we present the description of the chemical potential of 
QGP within the AdS/CFT correspondence. In Sect. 2.1.3 we remind the main facts 
about shock waves in AdSs related to the trapped surface formation. In Sect. 2. 1.4 
we describe in details the dual conjecture proposed in [p| . In Sect. 2. 2 we pay a 
special attention to the problem of regularization that appears within the wall shock 
waves approach. In Sect. 3 we present the phase diagram, temperature vs chemical 
potential, for QGP formed in the heavy-ions collisions by using the holographic 
approach with the central collision of charged shock waves. In Sect. 4 we present 
our calculations of the same problem by using the regularized version of the charged 
wall shock waves. We summarize our calculations in Sect. 5 and present here also 
further directions related to holographic description of quark-gluon plasma formed 
in heavy-ions collisions. 



2. Set up 

2.1 Dual Conjectures 

2.1.1 Black holes and AdS/CFT correspondence for strongly coupled QGP 

The idea of AdS/CFT applications to description of the QGP is based on the possibil- 
ity to make an one to one correspondence between phenomenological/termodynamical 
parameters of plasma - T, E, P, fj, - and parameters that characterize deformations 
of AdS 5 . In the dual gravity setting the source of temperature and entropy are 
attributed to the gravitational horizons. The relation between energy density and 
temperature typical for the BH in AdS according [^, is 

= STT^ 4 

16 G 5 y } 

In the phenomenological model of QGP, such as the Landau or Bjorken hydro- 
dynamical models |37], |3~8fl , the plasma is characterized by a space-time profile of the 
energy-momentum tensor T Atjy (x p ), fi,u,p = 0, ...3. This state has its counterpart on 
the gravity side as a modification of the geometry of the original AdSs metric. This 
follows the general AdS/CFT line: operators in the gauge theory correspond to fields 
in SUGRA. In the case of the energy-momentum tensor, the corresponding field is 
just the 5D metric. It is convenient to parameterize corresponding 5-dimensional 
geometry as 



^2 _ j2Qiiv(% p i z)dx^dx + dz ^ ^ 



z 2 
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that is the 5D Feffer man- Graham metric [39]. The flat case g^ v = r]^ parametrizes 
A0IS5 in Poincare coordinates. The conformal boundary of space-time is at z = and 

M<*) = ^ + A£iV) + --- (2-3) 

The AdS/CFT duality leads to the relation 

g^{x")-<T, u {x")> (2.4) 

Applications of AdS / CFT correspondence to hydrodynamical description of the 
GQP is based on the fact that the energy momentum tensor can be read off from the 
expansion of the BH in AdS$ metric ( |2.3| ) corresponding to simple hydrodynamical 
model 

< T, v >oc g$ = diag(3/4 1/4, 1/4, 1/4) (2.5) 



The BH in AdS§ in the Fefferman-Graham coordinates has the form (2^) with the 
following nonzero components of g fiu (x p , z) (see || and refs therein) 

C 1 -?) 2 ( z 4 \ 

9oo = - ) \{ , ^=1 + - (2.6) 
(l + |) V zo) 

A change of coordinates z = z/(l + z 4 /4) 1//2 transforms ( |2.2| ) to the standard metric 
form of the AdS-Schwarzschild static black hole 

z 2 ds 2 = -(l-^-] dt 2 + dx 2 + —^di 2 , (2.7) 
with z = z /\/2 being the location of the horizon. 

2.1.2 Chemical potential in QGP via AdS/CFT correspondence 

The Reissner-Nordstrom metric in AdS has the following form: 

ds 2 = -g{R)dT 2 + g(R)~ l dR 2 + R 2 dQ 2 D _ 2 , (2.8) 

where A is cosmological constant, A/3 = 1/a 2 , M and Q are related to the ADM 
(Arnowitt-Deser-Misner) mass m and the charge a 

M = M|m Q2= 4^. (210) 
6tt z 6 

a is a charge of the electromagnetic field (pure electric) with only one non-zero 
component 

A = A T dT = ( - + s)dr, (2.11) 
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\/3 Q 

here $ is a constant $ = where R+ is the largest real root of q(R). Ther- 

2 R\ K ' 

modynamics of the charged BH is described by the grand canonical potential (free 
energy) W = I/f3, the Hawking temperature T = 1/(3, and the chemical potential 
[ H| |4l| that are given by 

here -R+ is outer horizon, g(R+) = 0, J is given by the value of the action at (|2.9| ) 
and ( [2. 11| ). The relation with the first low of thermodynamics, d£ = TdS + fidQ is 



achieved under identification 

W = £-TS-$Q, £ = m, S = Q = q, // = $ (2.13) 

4G5 

Note that just the asymptotic value of a single gauge field component gives the 
chemical potential [p7| -pi 

H = lim A t (r) (2.14) 

r— >oo 

The QGP is characterized at least by two parameters: temperature and chemical 
potential. Generically speaking, quantum field theories may have non zero chemical 
potentials for any or all of their Noether charges. Within the AdS/CFT context two 
different types of chemical potential are considered, namely related to the R charge 
and to baryon number. 

Bar yon number charge can only occur when we have a theory containing funda- 
mental flavours. Introduction flavours into the gauge theory by means of a D7 brane 
leads to appearance of a U(Nf) global flavour symmetry. The flavour group contains 
a U(1)b, that is a baryon number symmetry, and for this baryon number one adds 
a chemical potential //& [p8fl . To calculate the free energy one has to calculate the 



DBI action for a D7 brane. Note that there is a divergence in formal definition, so 
we must go through the process of renormalization, see for example lectures and 
for yearly discussion (|3| . 



R charge chemical potential appears for SUSY models f34| . In the N = 1 case 
there is a U(l) R symmetry group. As to extended SUSY, say N=2, the quark mass 
term breaks R symmetry. 

The typical phase diagram the chemical potential vs the temperature is presented 



in Fig. D (the diagram is taken from [^7j). In the phase diagram: fj, q = j^-, jj q is the 



quark chemical potential and M oc m q is a mass scale defined as M = 2M q / vA and 
A = 9ym N c- 

2.1.3 Shock waves in AdS 5 

Shock waves propagating in AdS have the form 

, 9 r n—dudv + dx 2 \ + 6(x±, z)5(u)du 2 + dz 2 ,„*^s 
ds = L — \ ' y ' , (2.15) 
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B. 



Figure 1: Phase diagram from [27]: Quark chemical potential fi q /M q , in versus tempera- 
ture T/M. The red line separates the phase of Minkowski embeddings (small temperatures, 
small fj, q /M q ) from black hole embeddings (see details in |27]]). Figure (b) zooms in on the 
region near the end of this line. Different lines in B. correspond to different embedding 
geometries. 



where u and v are light-cone coordinates, and x± are coordinates transversal to the 
direction of motion of the shock wave and to ^-direction. This metric is sourced by 
the stress-energy momentum tensor Tmn with only one non-zero component 



rpSW 

uu 

and the Einstein E.O.M. reduces to 

3 



J uu (z,x±)S(u) 



where 



and 



Uh 3 ~ 7^) $ ( z > x -l) = -16nG s —J w (z,x±) 



$(z,x_l) = —</>(z,x±) 
z 3 d _ 1 d z 2 d 2 



(2.16) 

(2.17) 

(2.18) 
(2.19) 



Different forms of the shock waves correspond to different forms of the source 
J uu {z,x±). The most general 0(3) invariant shock wave in AdS located at u = 
corresponds to 

® m (z,x ± ) = F(q). (2.20) 

where q is the chordal distance 



(x\f + (x 2 ) 2 + (z- z ) 2 



Azz n 



In this case p, related to J uu as 



has the form 



-J uu (z,x ± ) = p(z,x±), 



p°( 3 \z,x ± ) = p(q), 



(2.21) 

(2.22) 
(2.23) 
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and the Einstein E.O.M takes the form 

(n Ha -^)F = -16icG s p(q) 



or explicitly 



q(q + l)F' q ' q + (3/2)(l + 2q)F' q - 3F = -16ttG 5 L 2 p{q), 



(2.24) 
(2.25) 



The shape of point shock wave F p is given by the solution to ( |2.17 ) with 

J uu = E5(u)5(z - L)5(x 1 )8(x 2 ) (2.26) 
and has the form is given by 



F p (z,x ± ) 



8L 2 G^Ez 3 



(2.27) 



(4 + (z- L) 2 ) 3 

This point shock wave shape is in fact equal to F p (q), $> pmnt (z,x±) = F p (q), that is 
a solution to ( 2.25 ) with 

E 5(q) 



L 3 y/q(l + q)' 



It has the form 



pp 



2G 5 E ( {8q 2 + 8q + 1) - A(2q + 1)^(1 + q) 



(2.28) 



(2.29) 



+ q) 

The shape of the charged point shock wave is a sum of two components 

F = F P + F Q (2.30) 
where F p is given by ( |2T29|) and F Q is the solution to ( ggg ) with 

5Q 2 1 5Q 2 1 



32 • 64 L 5 G 5 [q(q + l)] 5 / 2 vr24 • 64 L 5 + l)] 5 /2 ' 



(2.31) 



explicitly 



fQ _ 5G 5 Q 2 n lUq 2 + 16g - 1 + 128g 4 + 256g 3 - 64(2g + + l)ygj + g 
48^ 3 5(1 + 9)^(1 + 9) 

(2.32) 



To establish the connection with [22] let us note the relations of notations 

AG 5 E 



M 

Q 2 - 



3vr 
3tt 



(2.33) 
(2.34) 



and 



3ttM 2G 5 E 
2 a I notations ofQ ~~ 7^ I notations of IT and here 

57TQ 2 5G 5 Q 2 
g^ Q 3 I notations of [f22] 48L 3 'notations here 

More complicated shock waves in AdS and dS have been considered in [45 



(2.35) 
(2.36) 
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2.1.4 GYP Dual Conjecture 

Gubser, Yarom and Pufu (GYP) fllj have proposed the following dual to QCD 
holographic picture for colliding nuclei: 

• the bulk dual of the boundary nuclei is the shock waves propagating in AdS of 



the form ( 2.15 ) 



the bulk dual of two colliding nuclei in the bulk is the line element for two 
identical shock waves propagating towards one another in AdS 



2 2 — dudv + dx\ + z)8{u)du 2 + +4>2(x±, z)8{v )dv 2 + dz 2 . . 
ds — Li \ yZ.oij 

z 



when the shock waves collide in the bulk, a black hole should form, signifying 
the formation of a quark-gluon-plasma. 



To estimate the BH formation one usually use the TS technic |5(| ^TJ 2 . A 
trapped surface is a surface whose null normals all propagate inward |k|. There is 
no rigorous proof that the TS formation in asymptotically AdS spacetime provides 
the BH formation, however there is an common belief that trapped surfaces must lie 
behind an event horizon and that a lower bound on entropy Sa<is of the black hole 
is given by area of the TS, A tr a PP ed, 

O \ C — ^trapped / r> qo\ 

£>AdS ^ ^trapped = — TT^ [A. 66) 

To make the proposed duality prescription more precise one has to fix the rela- 
tions between the bulk parameters, G§,L,E and the phenomenological parameters 



of QGP. According to pT |, one of these relations is 

L 3 16E-T 4 11-16 



, , - 1-9 (2-39) 



The arguments supporting Q2.39| ) are following. Lattice calculations in QGP [pT 
have shown that ET A is a slowly varied quantity and 

ET 4 w 11. (2.40) 

Just to match the black hole equation of state ( 27T[) with ( 2.40 ), Gubser, Yarom 



and Pufu [11 1 have assumed ( 2.39 ). It is important to note that here is assumed an 
identification of the total energy of each nucleus with the energy of the corresponding 
shock wave. One can modify this identification and assume that only a part of energy 
of the gravitational shock wave is related to the total energy of nucleus. 



2 This estimation can be also performed using so-called capture arguments p2|, 
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To fix the dimensioless parameter EL the AdS/CFT dual relation (2.4) between 
the expectation value of the gauge theory stress tensor and the AdS$ metric defor- 
mation by the shock wave has been used [ITT] , 



(Tun) = t^t" I™ ^(z, x ± )5(u) (2.41) 

For the point shock wave $P° m * given by (|2.27|) , one gets the following stress tensor 
in the boundary field theory: 

< T -> = itp + ff + WW S(U) (2A2) 

The RHS of (|2.42|) depends on the total energy E and L, and L has a meaning of the 
root-mean-square radius of the transverse energy distribution. It has been assumed 



lT|] that L is equal to the root-mean-square transverse radius of the nucleons, that 



is in according with a Woods-Saxon profile for the nuclear density p4| , [55j is of order 
of few fm. In particular for Au it is equal to L « 4.3 fm. For Pb it is L w 4.4 fm. 

The RHIC collides Au nuclei, (A=197), at y/sj^ = 200 GeV. This means that 
each nucleus has energy E = 100 GeV per nucleon, for a total of about E = E be£ini = 
19.7 TeV for each nucleus. 

LHC will collide Pb nuclei, (A=208) at ^snn = 5.5 TeV, that means E = 

£beam = 570TeV. 



Estimations of [[TjJ for dimensionless values EL for Au-Au and Pb-Pb collisions 

are 

£L|A«_Au,^]7=200GeV ~ 43 X 1()5 » ( 2 ' 43 ) 

EL\ Au _ Au>VJ ^ =5 . 5TeY w 1.27 x 10 7 , (2.44) 

Note, that in [[TJ| has been proposed to tune the scale L or zq of the bulk colliding 
object to the size of the nucleus, or to the "saturation scale" Q s in the "color glass" 
models. 



Calculations in [TlJ show that in the limit of a very large collision energy E the 



entropy grows as 

£2/3 

3\ 1/3 



/ L 3 \ 1 

tapped ~ 7T {2EL) 2 ^ } (2.45) 

Considerations of off-center collisions of gravitational shock waves in AdS do not 
change the scaling E 2 ^ 3 . However, a critical impact parameter, beyond which the 
trapped surface does not exist has been observed [lj| (compare with result of p0|). 
Experimental indications for similar critical impact parameter in real collisions have 
been noted [|i~8| . 

The relation of the total multiplicity, Sqgp, given by experimental data, and 
the entropy produced in the gravitational waves collision in AdS^,, SAds has some 
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subtleties [21]. Phenomenological considerations |56, 



11], estimate the total mul- 
tiplicity Sqgp by the the number of charged particles N c h times the factor ~ 7.5. 



S, 



QGP 



7.57V, 



charged • 



(2.46) 



The trapped surface analysis does not give the produced entropy but it provides 
a lower bound 

Strapped < SacIS- (2.47) 



Taking into account that in calculations [|11|] the energy of the gravitational shock 
wave is identified with the energy of colliding ions and L with the nucleus size, one 
can introduce proportionality constants between these quantities to get 



■ St rCL pp e d *^ -^charged 



(2.48) 



where all proportionality factors are included into the overall factor M.. One can 

1/3 

take M. to fit the experimental data at some point. But the scaling Strapped <x s nn 



implied by ( [2.45 ) differs from the observed scaling, which is closer to the dependence 



S oc Sffx, that predicted by the Landau model |37|, see Fig.0. It is obvious, that 



if E < E max one can avoid a conflict between |TT| an d experiment, but if E can be 
arbitrary large the conflict takes place. 

In figure [| we plot the dependence of the entropy bound ( [2 .45 ) on the energy, 
together with the curve that schematically represents the realistic curve that fits 
experimental data |59|. We can see that by changing the coefficient A4 one can avoid 
the conflict only for energy up to some E max . The overall coefficient of the numerical 





5000 6000 



Figure 2: (color on-line) Plots of the total number of charged particles vs. energy. The 
red lines present the estimation (2.45). Plots A and B differ by the overall factor Ai. The 
blue lines correspond to the prediction of the Landau model and the dotted green lines 
schematically present the curves that fit experimental data. The dashed lines correspond 
to corrections to the GYP multiplicity via non-zero chemical potential, see Sect. 3. 



plot has been chosen in order to fit the RHIC data |59| . Their are indicated by dots 
in FigfJ. 



In the above estimation energy of each shock wave is identified with the energy 



of colliding beams. As has been noted in |L7j one can improve fit to the data by iden- 
tifying the energy of each shock wave with the fraction of the energy of the nucleus 
carried by nucleus that participate in the collision. This give an extra parameter to 
fit data. But still a conflict will arise at large energies. In paper [I7|| it has been 
proposed to cure the problem by removing a UV part of AdS bulk. In [ZS| shock 
waves corresponding to the BH with non-zero dilaton field [F24] were considered and 



it has been shown that lower bound on A/charged scales is rather closer to s 



1/4 
NN- 



2.2 Remarks about the regularization of TS calculations in the case of 
wall-on-wall collisions 



In [[Tj| has been proposed a much simpler dual description of the colliding nuclei that 
uses a wall-on-wall collision in the bulk. The Einstein equation for the profile of the 
wall shock wave [fTBI has the form: 



(d 2 z - -d z )4>{z) = J, 



WP 



J, 



WP 



z ) 



(2.49) 



To find a trapped surface that can be formed in the collision of two wall shock waves 
one needs to find a solution to the Einstein eq.( |2.49| ) that satisfies two conditions. It 
is convenient to write these conditions in terms of function ip(z) related to <fi via 



They have the form 



^\z a )j- 



2. 



z , 



0, 



(2.50) 

(2.51) 
(2.52) 



where z a , zt are supposed to be the boundaries of the trapped surface pq| . But as we 
will see in the moment, strictly speaking, one cannot call the solution to the equation 
( |2.49p with b.c. ( |2.51| ) and ( |2.52| ) the trapped surface, since by definition this surface 
supposed to be smooth and compact meanwhile the solution |TB| is non-smooth and 
noncompact. 

By this reason we call the solution found in a quasi-trapped surface. Let us 
remind the construction presented in 



In [ 50 1 , the solution to the Einstein equation Q2.49 ) is written in such a way that 
the property Q2.51 ) is satisfied automatically. This solution has the form 



Hz) 

^a(z) = 



-- ip a (z)<d(z - z) + ip b (z)Q(z - Zq) 
AGttE 



(2.53) 



zV*' 1 ZbZa 



L 4 (z b * - zj 



-12- 



Let us first note that solution (|2.53|) is not smooth. There is a non-smooth part 
of the solution ( |2.53| ) 



" K ( ~% ( T - % (Ta)) , (2-54) 



z \ z- zi 



Ti = z 4 Q(z -z) + zfo(z - zb) (2.55) 
T 2 = 4e(z - z) + z 4 Q(z - z ) (2.56) 



where 



Thus, in order to smooth the solution we have to smooth the function H. We can do 
it by performing the regularization of the Heaviside step function 

rv A ~ r arctan(#(2o-2;)) 3 1 

e(z - *) ~ "i = h - (2.58) 

7T Z 

0(* - z ) « T 2 = arCtan(i?(Z - Z ° ))3 + i : (2.59) 

7T 2 

and considering the regularized functions Ti and T 2 

~ 4 / arctan( J R(z - ^)) 3 A 4 / arctan (R (z - ^ )) 3 A 

^ + 2) +Z °[ n + 2) (2 - 60) 

~ 4 / arctan( J R(z - ^)) 3 1 \ 4 / arctan (R (z - z )) 3 A 

T 2 = ^ + -j+z ^ + -j (2.61) 

For derivatives we have 

^4, 3 9(, -.), ^1^4z3( arcta n(^ -,))3 + 7 r) ; 

az az 7r 

^4z 3 9(z-z ), ^ 4z3 (arc tan(^-z ))3 + 7r )^ 

az az 7r 

In Fig.|3| we present the derivatives of functions Tj, T 2 as well as derivatives of 
the smoothed functions T 1; T 2 . 



For R = 10 4 (see below) the differences between derivatives — - and their ap- 

dz 

proximations given by (|2.62|) and Q2.62j ) 

AiW = ^ - , A 2 (z) - — I — I (2X)\ ) 




AUz) = -A 2 (z) = -3- K — ^ h 3 - — ^ i ^— (2.65) 

KJ V> (1 + R^z -zf)-K (l + RS(z-z ) 6 )n K ' 

are of order ^ 10 3 fm 3 only in the interval z G [z' , Zq], Zq = 4.293 fm, zl = 4.307 fm. 

Indeed, in our consideration (spread case) the largest value of z a is 4.260706906 
fm and the smallest value of z b is 4.340400579 fm. At the points z' Q = 4.260706906 fm, 
4 = 4.340400579 fm the quantity Ai is less then < 5 ■ 10~ 6 fm 3 . 

At the points z' = 0.6948439783 fm, 4 = 1018.393720 fm the quantity Ai is 
less then < 2 • 10" 12 fm 3 . 




Figure 3: A. The functions — — (red line), — — 

dz dz 




B 



(blue line) . B. The functions 



dT 2 

dz 



appr 



(red line), 



dTi 



dz 



(blue line). The regularization parameter R = 10 at A and B cases. 



appr 



C. Functions — — (red line), — — 

dz dz 

parameter R = 10 4 . 



(blue line) and — — (green line) at the regularization 



appr 



The schematic picture of locations of roots and a region there |Aj(z)| > 10 -3 are 
presented in Fig.|j. We see that the difference A$ is not essential in location of the 
roots and we can use the approximations ( |2.62| ) and ( j2.63| ). 

The regularized version of the the function ip is 

Aeg = ^a{z)T 1 + ^) b {z)Y 2 . (2.66) 
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InE 



InE, . 



InE 



InE 




Figure 4: (color on-line) The schematic plots of locations of roots (solid black lines) 
dependent on the energy (in the logarithmic scale) and the location of differences 



( — ) 

\ J Jo 



> 10 



-3 



dz 

1 , 2 (the magenta shaded region) . The magenta solid line shows 



dz Japprl ~ 

the location of the wall. The dotted blue lines show location of zeros for the charged wall. 



Now one has to put conditions ( 2.52 ) on the regularized functions 



Zg d_ 

2L dz 

Zb d 

2L dz 



Ipreg 



Ipreg 



(2.67) 
(2.68) 



z=z b 



and find z a and z a from these conditions. However it is difficult to perform these 
calculations. Instead of finding z a from condition (|4.35|) we propose to use such 
regularization that does not change z a found from formal conditions (|2.52|). We 



can check that the formal z a in fact solves also the regularized condition if the 
regularization is smooth enough. So, we take z a and substitute it in the LHS of 
regularized condition (|4.35|) . We define 



F a . 



reg 



2L V dz dz 



b,reg 



Zb 

2L 



d^a r . d^ b 
in — ^—l 2 



dz 



dz 



= 2i, 



1+^1, 



-1 + S 2 . 



= 26 



We can calculate F a ^ reg . The deviation of F a ^ reg from 1 will show how the regu- 



larization changes conditions fl2.52|) . In the following table we present calculations 
of F a reg for the wide range of parameter of the theory. 

We choose the parameter R as minimally needed to make S\ and 82 negligible 
at energies 10~ 4 < E < 10 2 TeV. Using the direct numerical calculations we choose 
R = 10 4 . We perform numerical calculations at R = 10 4 and get the following table: 
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E, TeV 


Q, fm 1 / 2 


Za, f m 


Zb, fm, 


F a 


F b 


118.2 





0.04399350434 


4.015208900 • 10 6 


1.00000 


-1.00000 


30 





0.06948439782 


1.019088495 • 10 6 


1.00000 


-1.00000 


0.03 





0.6948439783 


1018.393720 


1.00000 


-1.00000 


0.00025 





4.260706906 


4.340400579 


0.99999 


-0.99999 



Thus, from the table evidently F a « 1, Ff, m — 1. 

As has been mentioned above, strictly speaking one may not consider infinite 
surface as a trapped surface of any kind. Nevertheless it is possible to assume that 
transversal size of colliding objects is finite but very large, and therefore boundary 
conditions do not affect the process of gravitational interactions of inner parts of 
sources. If we are interested only in the specific area of the formed trapped surface 
in respect to the unit of shock wave area, we may define it as 



A ~ lim 



L-»oo A source (L) 



(2.69) 



and the approximate equality takes place due to negligibility of boundary effects. As 
often happens, we can get answers for finite physical systems performing calculations 
for infinite non-physical objects. 



3. Holographic QGP phase diagram for the central heavy-ions 
collisions 

In this section we construct the phase diagram for TS formed in the central collision 
of two identical point-like charged shock waves ||22| . 

The profile of point-like charged shock waves in AdS is given by (|2.30| ) with 
( |2.29 ) and ( p.32j ). Existence of the trapped surface in the central collision of two 



point-like charged shock waves means the existence of a real solution, qo, to the 



following equation (see |22[ for details) 



2 2L 

i + v<Mi + qo) 
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The left hand side of (13.11) can be written as 



u s Af(L,M,Q 2 ,q) 



F(L, E, Q\ q) = v ' (3.2) 

V(a,q) 

The numerator M{L, E, Q 2 , q) contains just one term with dependence on Q 2 . This 
dependence is linear with a positive coefficient 

Af{a, M, Q 2 , q) = Af(a, M, q) + 15^ Q 2 . (3.3) 



The denominator in ( |3.2| ) does not take infinite values. To find solutions to ( |3.1| ) for 
the shape function given by ( |2.30| ) we can draw the function 

-Af(a, M, q) = -(512a V + 1280a V - 96Mvrag 2 + 1024a 3 g 3 - 96M7iaq 

+ 256a 3 g 2 ), (3.4) 

and see where this function can be equal to a given value 15Q 2 ^. 

In order to find the maximal allowed Q 2 at which solution to (|3.1|) still exists we 
find the maximum of function M for fixed energy, 

dAf(a,M,q) l 

^ \q=qmax — U \°-°) 

and the value 

^-Af(a,M,q)\ q=qmax 



defines Q 2 max . 

Let us remind that we are working in physical units and we use the following 
notations (KM) and ( Egg) : M = and Q 2 = 5 ^ n . 

37T 37T 

Results of calculations are presented in Fig. H. 



-17- 




20 40 60 80 100 



E, Tev 



Figure 5: The allowed zone for the trapped surface formation is under the line on the 



diagram. The plot has been constructed by using formulas from 22 . 



To estimate corrections to GYP multiplicity due to non-zero chemical potential, 
we use formula (3.17) from [22]. In notations admitted in this paper, ( p.29| ) and 
([2.32 ), the formula has the form 



.4 



AdS 5 



In Fig|| we show the entropy, A^dSs, f° r Qn = and Q n ^ 0. The blue line 
represents Q n = 0. The red line represents Q n = 2 • 10 6 . We see that the deviation 
form the GYP multiplicity is essential for small energies and is almost neglectful for 
large energies. 




E, TeV 



Figure 6: The function A^ds^i a t Qn = 0, (blue line) and Q n = 2 • 10 6 /m 1//2 (red line). 
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4. Holographic QGP phase diagram in the wall- wall dual model 
of heavy- ions collisions 



4.1 Charged wall as a dual model for a heavy-ion with non-zero chemical 
potential 

Let us note that the form of the J^ p in fl2.49|) can be obtained by spreading out the 
energy-momentum tensor of an ultrarelativistic point, i.e J uu in the form (|2.22| ) with 
p(q) given by eq.( |2.28|) , over the transversal surface. 

The Einstein equation for the charged wall (membrane) has the form 

(d 2 z - \d s )4>(z) = -16vrG 5 {JZ P + JZ Q (Q, z)) . (4.1) 

where is given by ( |2.49| ) and we suppose that JZ Q (.Qi z ) can De obtained in 
the similar way by spreading the energy-momentum tenzor of the ultrarelativistic 
charged point over the transversal surface. In the previous calculations: 



jWQ _ Im Juu^X^ ^ 



here the subscript "pQ" means the electromagnetic part of the energy momentum 
tensor of the charged point particle and "Vx±" means that we integrate over the 
induced metrics on the orthogonal surface Ai. 
For this purpose we take 

J%(z,zo) = ~/ Q (4-3) 

where ff® is given by ( |2.31|) , and according to our prescription fl4.2|) we integrate 
over all transversal coordinates 



The result is 



J* = XJ (4.5) 



where 



3 _ ^ ^ _ 4-3^-yg + ^ ( , 6) 
x _ _^_Ql_ _ _^_Q^_ ( 47 \ 

256 ttL 6 256 L 6 1 ' ; 

We see divergency at z = zq, as it should be for the energy- momentum ten- 
sor of a charged plane. We introduce regularization by adding the e factor in the 
denominator. 
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4.2 Charged wall-on-wall collision as a dual model for heavy-ions collisions 
with non-zero chemical potential 

To find the TS formation condition in the wall-wall collision one has to solve Einstein 
equation 

{dl - - z d z )4>{z) = -16nG 5 (Jg(z) + Jl W (Q,z)) , (4.8) 

(4.9) 



E ~ 3 



J uu (z) = ^S(z-z ), (4.10) 

lQW(n 128* z A (-z 2 + 3 z 2 ) 6(z -z)+ 4 {-3z 2 + z 2 ) 6{z - z ) 

J uu iQ, z ) ~ -^—zzo : — — 2 2 , 3 

3 {-z 2 + Zq + e 2 ) 

(4.11) 

with the following boundary conditions 

1) <f>(z a ) = <f>{z b ) = 0, <f) a {z ) = Mzo) (4.12) 

2) (V(^)f)=2, (>(* 6 )|)=-2, (4.13) 

where z a and z b are the boundaries of the TS and if) is related to 

<f>{z) = ^. (4.14) 

We search for a solution to the Einstein equation with a charged source in the 
form of the sum of the "neutral" solution and a correction proportional to Q 2 

<P = <Pn + (f> q (4.15) 

here n denotes the solution of the neutral case. 

As in the neutral case it is convenient to consider domains z < z , z > z 
separately 

^qzo>zi Zq > Z, 
fiqzyzo-, Z > Zq 



(4.16) 



and we have 



(& i« n v 128 ^{-z 2 + ?,zq 2 ) 

{dl - -d z )(f) qzo>z = -16ttG 5 X—zzo- — z > z; (4.17) 

z 6 {-z 2 + Zq* + e 2 ) 

{dl - -d z )<P qz >z = -16nG s X^zzo f°[~ 3z2 2 +4 L ,*>«>■ (4-18) 
z 3 {—z 2 + Zq — e 2 ) 



Solutions to (|4.17|) and ( |4.18|) can be presented as 

C 2 NLz z 3 

z 4 {-z 2 + z 2 + e 



1pqz >z = Z 3 C X + — — — — - , Zq > Z, (4.19) 



5 



%z> zo = — + z 3 C 4 + 2V z>z (4.20) 

Z 4:Z { — Z 2 + Zq 2 + e 2 ) 
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Here N = f ^ff- The first two terms in (|4l9|) and ( ^20|) are solution to the Lin 
and Shuryak equation (55) in fl8| . If one assumes that they satisfy condition 1, i.e. 

*Pn(Za) = tpn{Zb) = 0, Vw(^o) = VVlft^o), One gets |60| \ 



1 / 



C Kr-- 



C 



4vrG,£ \ z\ 



L 4 



^ < -20 



^6 



raft 



4ttG s £ 



-1U 



L 4 



2 < 2 



3~3 



2« 2, 



(4.21) 



In the neutral case one find z a and z b from the 2-nd condition (ip' n a( z a)^) 



(Vk(*)f 



-2, here z a and z& are the boundaries of the TS. 



As to ( g[Ig ) and ( p~20| ), choosing 



C 3 



jVLz 
4(^-z 2 )' 
NLzq 



Co. = 0, 



c 4 = o, 



(4.22) 
(4.23) 



we obtain 



NLzqZ 3 



-4 + z 2 -e 2 



4 (-^ + z 2 + e 2 )( _ z 2 + z 2 



2V 



6g 



NLzq 



-zl + z 2 + e 2 



Z < Z 



Z < Z 



(4.24) 



4z (-^+^- e 2) ( _^ + ^)' 

Note that for the constructed solution the condition ip{z a ) = ip(zb) = is satisfied 
automatically. 

The second requirement (|4.13[) gives 



8nG 5 E(4-4)z* a N z,z\ 
8vrG 5 E (4 - 4) zl N 4z h 

+ t / .2 , = _1 ; ( 4 - 26 ) 



L5(4-z 4 ) 4(-^ 2 + z 2 )2 

These equations do not have analytical solutions and we treat them numerically. 

Roots of system (|4.25|) , (fL~2"ED could not be found analytically since these equa- 
tions are equivalent to polynomial equations on z a and zj, of a high degree (> 4). So 
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we take zq = L and analyze the following system numerically 



F a = 



F h EE 



8kG 6 E(4-4)z* IOttGsQ 2 



Z Z a 



z o( z b 



4) 4 



z o( z b 



+ 



10 ttG 5 Q 2 



-4 + 4? 

4 z b 



[-4 + 4) 2 



(4.27) 
(4.28) 



To show the movement of roots of equations (|4.27|) and (|4.28|) we suppose that 
z for given Q is already known and represent function F a (z a , z ) as function of z a in 
Fig. [7]. In the similar way, supposing that z a is already known we represent function 
Fb(z a , z ) as function of Zb in Fig. E. 



50? ■ 10 :i 




Fa 
1.2 




20V*f ■ 10 3 



400 ■ 10 3 



0.05 z a B. 

Figure 7: The plot of F a (z a , Zb) as a function of z a for fixed Zb near the root z a = z a (E) 
at E = 118.2 TeV. Figure (B) zooms in the region of small F a and small z a . 




Figure 8: The plot of Fb(z a ,Zb) as a function of Zb for fixed z a near the root Zb(E) at 
E = 118.2 TeV. Figure B zooms in the region of small negative Fb and presents Zb in the 
logarithmic scale. 
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In Fig|| we show the charge flows of the roots. Different lines correspond to 
different energies. We see that the flows go to zq and reach the line z = zo for 
Q — Qcr- in FigJIDIwe draw the corresponding flow for physical parameters. 




Figure 9: The schematic picture of charge flows The magenta solid line shows the position 
of the wall. We see that the positions of points z a (Q) and z b {Q) move to the point z = zq 
when we increase Q. For Q — > Q cr (E) the segment [z a (Q), Zb(Q)] shrinks to zero. 



1, x 10 6 
9, x 10 5 
8, x 10 5 
7, x 10 5 
6,x 10 5 
5,x 10 5 
4,x 10 5 
3,x 10 5 

2, x 10 5 
Lx 10 5 



Figure 10: A. The charge flows of the root z a (Q) for E = 1.97 TeV, 3TeV, 9TeV, and 
24TeV. B. The charge flows of the roots z b {Q) for E = 1.97 TeV, 3 TeV. 

4.3 Comparison of results. 

It is interesting to compare the phase diagrams, the energy (temperature) E vs the 
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charge (chemical potential) Q, corresponding to the pointlike charge and the spread 
charge. Results of these calculations are collected in the table below and presented 
in Fig.O. We see that this two phase diagram are almost the same. 



E (TeV) 


118.2 


60 


30 


6 


3 


0.6 


0.06 


0.03 


0.0003 


0.00025 


Qcr, point 
Qcr, wall 


25649.6 
47500 


14577.2 
27000 


8180.6 
15170 


2138.7 
3950 


1199.9 
2220 


313.3 
570 


45.6 
80 


25.4 
40 


0.43 
0.15 


0.37 




HQft) Q„,/m 1/!1 




Figure 11: A. The phase diagram the logarithm of Q n vs the logarithm of E at large E. 
B. The phase diagram E vs Q n for small E and small Q n . The blue lines correspond to 
the pointlike charge and the red lines to the spread charge. The zones above the lines are 
forbidden for black holes production for corresponding E and Q. 



From Fig. |TT] it is evident that the two lines, the red and the blue ones, have a 
cross point. We represent the cross point in natural and logarithmic scales in Fig. 
12. 



HQ,,) Qn./m 1 ' 2 




Figure 12: The cross point of two diagrams in logarithmic and natural scales. 



4.4 The square trapped surface calculation 



Following [|18[ we calculate entropy lower bound as "the area of the trapped surface" 
per an unite square of the wall 3 using the formula: 



S 



2A J ^dzd 2 x± 
S 

Jd 2 x ± ~ 4G 5 



2G 5 
L 3 (\ 



(4.29) 
(4.30) 



In the absence of transverse dependence one ignores x\ in (|4.29| ). QO0| ) measures 
entropy per transverse area. 

The trapped surface decreases with growth of a charge. The corresponding 
graphical representations are in Fig. |13|. 




3000 Q„/0F 




30000 40000 Q n j^TT 



Figure 13: The dynamics of the trapped surface area s(<5n/v / ^') at .E7 = 6TeV, E 
118.2 TeV. 



In Fig.[I4] we show the entropy per volume given by ( [4.30| ) as function of energy 
for different Q. This plot is similar to the plot presented in Fig. [5| We see that the 
influence of the chemical potential on the multiplicity is essential for small energies 
and is almost neglectful for large energies. 



4.5 Remarks about the regularization 

The regularized version of the the function ip is 

where ip a {z) and ipb{z) define the function ip without regularization, 

ip = tp a (z)O(z -z) + i) h {z)Q{z - z ) 



(4.31) 



(4.32) 



3 We put "area" and "trapped surface" in quotation marks since in the strict notions of the 
trapped surface it has to be smooth and compact. In our case it is not smooth and it does not have 
finite area, one can only assume this properties after regularization 
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s, fm 




Qn = 



Qn = 5000,/? 



}„ = 2000,/? 



E, TeV 



Figure 14: The red line corresponds to the case Q n = Ofm 1 / 2 , the blue to the case 
Q n = 2000^ fm 1 / 2 , the black to the case Q n = 5OOO0F fm 1 / 2 . 



and 



^E[^-l)z h W[p- z ^) io Q 2 Gnz zH-z a 2 + z 



L 4 (z b * - z a *) 



4GnEl3$-l)z a W ,JL 



L* (z b * - z a *) 



+ 



3 L 5 (-z 2 + zo 2 ) {-z 2 + z 2 ) 

10 Q 2 GlT Zp 5 (~Z b 2 + Z 2 ) 

3 Jj> z (_ Z 2 + ^2) (_^2 + ^2) 



arctan (R (zq — z)) 1 
1 1 = Z h o 

arctan [R (z — z )) 1 

F 2 = Z h 9 

7T Z 



Now one has to put conditions ( |2.52| ) on the regularized functions 

z a d 



2L dz 



Ipreg 



(4.33) 
(4.34) 

(4.35) 



However it is difficult to find z a from the condition ( 4.35|) . fnstead of finding z a from 
the condition (|4.35|) we propose to use such regularization that does not change z a 
found from the formal conditions ( [2.52 ). We can check that the formal z a in fact 
solves also the regularized condition if the regularization is smooth enough. So, the 
take z a and substitute it in the LHS of the regularized condition ( f4.35|) . We define 



F a . 



reg 



'dipa T dj) b ^ 
2L V dz dz 
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b,reg 



Zb_ ( <tya_ T 

2L \ dz 1 cfe 2 



-1. 



We can calculate F a ^ reg and F^ reg . In the following table we present calcula- 
tions of F a>reg and Fb jreg for the wide range of parameter of the theory. Results of 
calculations at R = 10 4 are presented in the following table: 



E, TeV 


Q = Qn/V^, M 1/2 


Za, fm 


Zbjm, 


F a 


F b 


118.2 


40000 


0.04928014740 


4.015208864- 10 6 


0.99997 


-1.00000 


3 


15000 


0.08847525298 


1.019088359- 10 6 


1.00000 


-1.00000 


0.03 


40 


0.7861838575 


1017.792389 


1.00000 


-1.00000 



Thus, from the table it is evident that F a pa 1, F b m — 1. 

5. Conclusion 

5.1 Summary 

In this paper we have constructed the phase diagram of the quark gluon plasma 
(QGP) formed at a very early stage just after the heavy ion collision. In this con- 
struction we have used a holographic dual model for the heavy ion collision. In this 
dual model colliding ions are described by the charged shock gravitational waves. 
Points on the phase diagram correspond to the QGP or hadronic matter with given 
temperatures and chemical potentials. The phase of QGP in dual terms is related 
to the case when the collision of shock waves leads to formation of trapped surface. 
Hadronic matter and other confined states correspond to the absence of trapped 
surface after collision. 

Multiplicity of the ion collision process has been estimated in the dual language 
as an area of the trapped surface. We have shown that a non-zero chemical potential 
reduces the multiplicity. To plot the phase diagram we use two different dual models 
of colliding ions. The first model uses the point shock waves and the second the wall 
shock waves. We have found qualitative agrement of the results. 

A special attention has been devoted to a regularization procedure for calcula- 
tions performed for wall shock waves. On the one hand technically these calculations 
are essential simpler, but on the other hand, this approach, strictly speaking, is incor- 
rect and requires a regularization. We have shown that a natural regularization does 
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exist. Moreover, the proposed regularization does not make calculations to be more 
complicated as compare with naive (direct) calculations. This open new possibility 
for simple calculations for wall shock waves bearing nontrivial matter charges. 

5.2 Further directions 

Head-on collisions of point charged shock waves have only two parameters. In the 
dual language they correspond to energy and chemical potential per nucleus. Off- 
center collisions are also specified by the impact parameter and the change of this 
parameter can be associated with a dual change from non-thermal peripheral to 



thermal central collisions [[L8|j . However, this is still an oversimplification of the 
problem. The physics of heavy-ion collision in RHIC is richer and as indicated in 
[p!8| , |60| , rapid equilibration and hydrodynamical behavior experimentally observed 
at RHIC for collisions of two heavy ions such as AuAu, does not have the place 
for deuteron-Au collisions at the same rapidity. Maybe it is too naive to believe 
that the simplest shock wave related by a boost to the Schwarzschild black hole 
in AdS can mimic the nuclear matter in the colliders. However this simple shock 
wave in fact reproduces the interaction of a relativistic quark with gravity and by 
this reason, may be considered as a simplest candidate to mimic the nuclear matter 
within the holographic conjecture. One can try to associate different nuclei with 
different forms of shock waves. Let us remind in this context that the form of the 
shock wave follows from the eikonal approximation of the gravity-quark interaction 
in 5-dimension |)T, 5B|. The presence of the electromagnetic field or other fields as 



well as any improvements of the eikonal approximation for sure changes the form of 
the shock waves and it would be interesting to see holographical consequences of this 
consideration. 

1/3 

The obtained lower bound on A^ har g ed scales as s^ N , which is a faster energy 
dependence than the s]^- scaling predicted by the Landau model |37j and largely 



obeyed by the data. If one has a priory a restriction on allowed energy then one 
can fit constants to guaranty that the experimental data are above the AdS bound. 
Note that taking into account the chemical potential permits to increase the allowed 
energy. However one cannot expect to much from the chemical potential corrections. 
The relevant chemical potential for baryon number is not expected to be large, i.e. 
fis ~ 30MeV or Hb/T ~ 0.15 for recent experiments at RHIC pSI and so any effects 
will be limited. However, as has been mentioned in the text, the relation between 
the value of chemical potential and the value of the 5-dimentional charge is in our 
disposal and we can assume a huge ratio of them. 

It would be also interesting to try to use plane gravitational waves in AdS 5 to 
describe nonperturbative stages in the gauge theories and collisions of these waves 
to describe the QGP formed in the heavy ions collisions. In the plane case, the 
Chandrasekhar-Ferrari-Xanthopoulos duality between colliding plane gravitational 
waves and the Kerr black hole solution, has been used as a model of the BH formation 
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j33[l . It would be interesting to generalize this duality to the AdS case. This may get 
a new insight to a possible dependence of multiplicities on the rapidity. 
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